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ABSTRACT  NOMENCLATURE 


The  wake  behind  a  two-bladed  model  rotor  in 
light  climb  was  measured  using  particle  image  ve- 
locimetry,  with  particular  emphasis  on  the  develop¬ 
ment  of  the  trailing  vortex  during  the  first  revolution 
of  the  rotor.  The  distribution  of  vorticity  was  dis¬ 
tinguished  from  the  slightly  elliptical  swirl  pattern. 
Peculiar  dynamics  within  the  “void”  region  may  ex¬ 
plain  why  the  peak  vorticity  appeared  to  shift  away 
from  the  center  as  the  vortex  aged,  suggesting  the 
onset  of  instability.  The  swirl  and  axial  velocities 
(which  reached  44  and  12  percent  of  the  rotor-tip 
speed,  respectively)  were  found  to  be  asymmetric  rel¬ 
ative  to  the  vortex  center.  In  particular,  the  axial 
flow  was  composed  of  two  concentrated  zones  mov¬ 
ing  in  opposite  directions.  The  radial  distribution 
of  the  circulation  rapidly  increased  in  magnitude  un¬ 
til  reaching  a  point  just  beyond  the  core  radius,  after 
which  the  rate  of  growth  decreased  significantly.  The 
core-radius  circulation  increased  slightly  with  wake 
age,  but  the  large-radius  circulation  appeared  to  re¬ 
main  relatively  constant.  The  radial  distributions 
of  swirl  velocity  and  vorticity  exhibit  self-similar  be¬ 
haviors,  especially  within  the  core.  The  diameter  of 
the  vortex  core  was  initially  about  10  percent  of  the 
rotor-blade  chord,  but  more  than  doubled  its  size 
after  one  revolution  of  the  rotor.  According  to  vor¬ 
tex  models  that  approximate  the  measured  data,  the 
core-radius  circulation  was  about  79  percent  of  the 
large-radius  circulation,  and  the  large-radius  circu¬ 
lation  was  about  67  percent  of  the  maximum  bound 
circulation  on  the  rotor  blade.  On  average,  about  53 
percent  of  the  maximum  bound  circulation  resides 
within  the  vortex  core  during  the  first  revolution  of 
the  rotor. 


‘Presented  at  the  American  Helicopter  Society  59th  An¬ 
nual  Forum,  Phoenix,  Arizona,  May  6-8,  2003.  Copyright 
©2003  by  the  American  Helicopter  Society  International,  Inc. 
All  rights  reserved. 


a  speed  of  sound,  meters  per  second  (m/s) 

A  area  enclosed  by  circulation  path,  m2 

c  chord  length,  m 

Cq  torque  coefficient,  torque/p7rfi2f?5 

Ct  thrust  coefficient,  thrust/p7rft2R4 

dh  distance  between  velocity  peaks,  horizontal  cut 

dv  distance  between  velocity  peaks,  vertical  cut 

M  hover- tip  Mach  number,  RQ./a 

N  number  of  points  to  be  averaged 

N),  number  of  rotor  blades 

R  radius  of  rotor-tip  path,  m 

Re  hover-tip  Reynolds  number,  cRQ/v 

r,  9 ,  z  coordinates  relative  to  rotor  axis 

u,v,w  velocity  components,  m/s 

uc,wc  vortex-convection  velocities,  m/s 

ur,wr  velocities  relative  to  uc,wc,  m/s 

a  angle  of  position  vector,  degrees  (deg) 

1 3  angle  of  u  +  w  vector,  deg 

T  circulation,  m2/s 

rc  circulation  at  core  radius,  m2/s 

F*  total  vortex  circulation  in  far  field,  m2/s 
6  turbulent/laminar  viscosity,  e/v  ' 

e  eddy  viscosity,  m2/s 

C  radius  measured  from  void  center,  mm 

Cc  radius  of  vortex  core,  millimeters  (mm) 

Cco  value  of  Cc  at  trailing  edge,  mm 

C  C/Cc 

9  local  angle  of  rotor  blade,  deg 

k  radius  measured  from  vorticity  centroid,  mm 

A  dimensionless  hover-inflow  velocity,  Wi/QR 

i/  kinematic  molecular  viscosity,  m2/s 

p  fluid  density,  kg/m3 

<t  rotor  solidity,  Nbc/rR 

swirl  velocity,  m/s 

ipc  swirl  velocity  at  core  radius,  m/s 

<P  V/Vc 

ip  azimuthal  wake  age  measured  from  blade,  deg 

vorticity  normal  to  r  —  z  plane,  per  second 
fi  rotor  revolutions  per  minute,  rpm 
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INTRODUCTION 

A  rotor  wake  follows  a  trajectory  that  repeat¬ 
edly  passes  near  its  own  generating  rotor  blades  be¬ 
fore  finally  being  convected  downstream.  As  a  con¬ 
sequence,  the  rotor  wake  not  only  has  a  considerable 
influence  on  the  performance  of  the  rotor,  but  it  is 
principally  responsible  for  the  high  levels  of  noise  and 
vibration  that  result  from  encounters  between  the  ro¬ 
tor  blades  and  their  collective  wake  (ref.  1).  Because 
of  the  complexity  of  the  rotor  wake,  especially  con¬ 
cerning  the  intensity  and  size  of  the  dominant  trail¬ 
ing  vortex  and  its  increasingly  irregular  trajectory, 
computational  modeling  of  the  flow  field  continues 
to  be  a  challenge.  Reliable  experimental  measure¬ 
ments  are  therefore  needed  to  support  such  efforts. 

As  new  experimental  techniques  emerge,  the  ro¬ 
tor  wake  is  frequently  revisited  in  order  to  extract 
more  detailed  information  about  the  structure  of  the 
trailing  vortex  (refs.  2  -  15).  For  example,  not 
only  does  the  stereo  particle  image  velocimetry  (PIV) 
measurement  technique  render  all  three  components 
of  velocity  over  a  fairly  large  cross  section  of  the  flow, 
but  the  vector  fields  from  a  sequence  of  images  taken 
during  successive  revolutions  of  the  rotor  can  also  be 
conditionally  ensemble  averaged  to  remove  the  ef¬ 
fects  of  vortex  wander  (ref.  14).  Given  this  opportu¬ 
nity  to  acquire  data  more  efficiently,  highly  resolved 
(in  space  and  time)  data  on  the  development  and  ag¬ 
ing  of  the  trailing  vortex  is  now  possible.  In  addition 
to  enhancing  our  physical  understanding  of  the  rotor 
wake,  specific  characteristics  (such  as  composition, 
size,  intensity,  and  location)  of  the  trailing  vortex  at 
different  wake  ages  will  enable  computational  fluid 
dynamics  (CFD)  models  to  be  validated  on  a  more 
detailed  level. 

Many  questions  have  been  raised  concerning  un¬ 
expected  features  appearing  within  the  trailing  vor¬ 
tex  (ref.  3),  the  composition  of  the  vortex  relative 
to  blade  loading  (ref.  2),  vortex  definitions  (ref.  16), 
opinions  dealing  with  the  development  and  demise  of 
vortex  systems  (ref.  17),  and  the  appropriateness  of 
simple  vortex  models  (ref.  18).  In  response  to  these 
questions,  an  experiment  was  performed  to  provide 
data  that  focuses  specifically  on  cross  sections  of  the 
trailing  vortex  during  a  complete  revolution  of  the 
rotor.  The  data  are  of  sufficient  quality  to  permit 
various  modeling  uncertainties  to  be  resolved.  This 


paper  describes  the  test  facility,  model,  flow  mea¬ 
surement  technique,  and  the  data-acquisition  and  - 
processing  procedure.  A  discussion  of  the  results  is 
provided  along  with  specific  findings.  A  considerably 
more  detailed  presentation  is  given  in  reference  19. 

TEST  DESCRIPTION 

Model  Rotor  Setup 

A  model  rotor  was  placed  in  the  center  of  the 
settling  chamber  of  the  7-  by  10-Foot  Subsonic  Wind 
Tunnel  at  Ames  Research  Center,  and  mounted  hor¬ 
izontally  with  its  axis  of  rotation  coincident  with  the 
centerline  of  the  test  section.  Aside  from  the  influ¬ 
ence  of  the  tunnel  walls  (approximately  4.6  m  away) 
and  the  contraction  into  the  test  section  (with  an 
area  ratio  of  14:1),  the  rotor  wake  was  allowed  to 
pass  through  the  test  section  without  interference. 
A  curtain  was  drawn  across  the  air-exchanger  sec¬ 
tion  to  prevent  the  wake  from  recirculating  around 
the  tunnel  while  fresh  stagnant  air  was  drawn  by  the 
rotor  into  the  tunnel  circuit  (as  required  by  mass 
continuity) . 

The  rotor  used  in  this  experiment  was  a  1/7- 
scale  model  of  the  full-scale  AH-l/OLS  rotor  (fig.  1). 
The  rotor  consisted  of  two  blades,  each  with  a  chord 
of  103.9  mm,  maximum  thickness  of  10.2  mm,  and  a 
span  of  957.6  mm  measured  from  the  axis  of  rotation 
to  the  blade  tip.  Beginning  at  25-percent  radius,  the 
twist  of  each  blade  was  — 0.106°/centimeter  (cm),  so 
that  the  total  variation  in  geometric  twist  along  the 
span  was  7.55°  (from  the  center  of  rotation  to  the 
rotor  tip).  Collective  pitch  settings  were  referenced 
to  the  geometric  pitch  angle  at  75-percent  radius. 
Therefore,  a  collective  pitch  setting  of  zero  meant 
that  the  blade  pitch  was  about  5.0°  at  the  root  (held 
fixed  for  locations  less  than  Rj 4)  and  —2.5°  at  the 
tip.  The  rotor  hub  had  a  built-in  blade  precone  of 
1.30°. 

Stereo  PIV  Setup 

The  primary  elements  of  the  PIV  system  con¬ 
sisted  of  a  pulsed  laser,  optics  to  spread  the  beam 
into  a  thin  sheet,  charged  coupled  device  (CCD) 
cameras  to  record  the  images,  and  a  calibration  tar¬ 
get.  The  plane  of  the  light  sheet  was  horizontal  and 
passed  just  below  the  trailing  edge  of  the  rotor  when 
the  rotor  was  at  zero  azimuth  (fig.  2).  In  order  for 
both  cameras  to  receive  close  to  the  same  amount 
of  forward-scattered  light,  two  opposing,  coplanar 
light  sheets  were  established  using  synchronized  120 
millijoule  (mJ)  Nd.YAG  lasers.  Cameras  and  lasers 
were  triggered  off  of  a  1/rev  reference  signal  that  was 
keyed  to  the  rotor.  This  signal  was  obtained  by  di¬ 
recting  a  small  continuous-wave  laser  beam  onto  a 


small  mirror  attached  to  the  hub  of  the  rotor  and 
intercepting  the  reflected  beam  with  a  fast-response 
photodetector.  Acquisition  software  was  used  to  de¬ 
lay  the  reference  signal  in  order  to  obtain  data  at 
different  wake  ages. 

Images  were  acquired  with  8-bit  CCD  cross- 
correlation  cameras,  each  having  a  sensor  array  of 
2K  x  2K  pixels  (nominal).  Each  pixel  measured  7.4 
microns  (/^m)  on  a  side.  The  cameras  were  operated 
in  a  double-exposure  mode  to  acquire  two  noninter¬ 
laced,  full-frame  images  during  a  single  frame  inter¬ 
val.  The  lenses  were  remotely  translated  to  focus  on 
the  centerline  of  the  image  area,  and  the  sensor  (lo¬ 
cated  inside  the  camera  body)  was  rotated  about  its 
centerline  to  satisfy  the  Scheimpflug  condition1. 

The  lasers,  sheet-forming  optics,  and  cameras 
were  rigidly  mounted  on  a  common  structure,  which 
could  be  traversed  in  three  orthogonal  directions. 
This  arrangement  allowed  the  cameras  to  be  focused 
only  on  the  trailing  vortex  region  of  the  flow,  and 
then  traversed  as  required  to  follow  the  position  of 
the  vortex  at  different  wake  ages. 

Proper  seeding  of  the  flow  is  critical  for  obtain¬ 
ing  accurate  PIV  measurements.  The  seed  particles 
must  be  evenly  distributed  and  of  sufficient  density 
to  define  the  flow  without  altering  its  physical  prop¬ 
erties.  The  particles  must  also  be  small  enough  to 
accurately  follow  the  flow  (especially  challenging  in 
accelerating  flows),  yet  large  enough  so  that  they 
scatter  a  sufficient  amount  of  light  to  be  detected. 
The  particle  generator  used  in  this  test  employed 
an  inert  gas  to  atomize  a  nontoxic,  pharmaceutical- 
grade  mineral  oil.  The  mist  was  vaporized  and  then 
condensed  before  being  released  into  the  flow.  The 
particle  size  was  estimated  by  the  manufacturer  to 
be  less  than  0.5  pm.  The  particle  generator  was  lo¬ 
cated  upstream  of  the  settling  chamber  so  that  the 
particles  would  mix  with  the  air  entering  the  tunnel 
circuit  before  being  drawn  into  the  wake  of  the  rotor. 

PROCESSING  WAKE  IMAGES 

To  obtain  good  resolution  of  the  trailing  vortex, 
a  calculation  grid  containing  92  horizontal  nodes  by 
79  vertical  nodes  was  constructed  over  the  region  of 
significant  interest  in  the  flow  field.  The  area  cov¬ 
ered  by  the  grid  resulted  in  approximately  12  pixels 
between  nodes  in  both  directions.  The  interroga¬ 
tion  window  was  set  at  22  pixels  on  a  side,  giving  a 
45-percent  overlap.  The  physical  area  measured  in 
each  image  zone  was  168  mm  (with  an  axial  incre¬ 
ment  of  1.8  mm)  by  106  mm  (with  a  radial  increment 
of  1.4  mm).  The  1.8-  by  1.4-mm  increments,  which 
defined  the  calculation  step  sizes  in  the  axial  and 

Uniform  focus  requires  that  the  light  sheet,  camera  lens, 
and  camera  sensor  planes  intersect  along  a  common  line. 


radial  directions,  also  defined  the  resolution  of  the 
measurements.  Therefore,  the  interrogation  window 
(containing  484  pixels)  covered  an  area  measuring 
2.5  mm  in  the  radial  direction  by  3.4  mm  in  the  ax¬ 
ial  direction. 

Consider  an  interrogation  area  in  which  10  par¬ 
ticle  pairs  (the  minimum  number  required  to  yield  a 
valid  resultant  vector)  have  been  registered  on  a  22- 
by  22-pixel  image  of  a  small  region  in  a  swirling  flow. 
The  software  derives  a  representative  in-plane  veloc¬ 
ity  vector  based  on  a  second-order  calculation  of  all 
the  individual  particle  displacements  within  that  in¬ 
terrogation  area.  Note  that  the  effective  “probe  size” 
is,  therefore,  less  than  the  interrogation  area.  All  the 
rotational  velocity  components  within  the  vortex  are 
accounted  for,  but  they  are  represented  in  a  collec¬ 
tive  sense  over  contiguous  regions  defined  by  the  size 
of  the  interrogation  area.  For  an  initial  core  diameter 
of  10  mm,  the  core  would  have  been  covered  by  about 
nine  interrogation  areas  (78.5  mm2/8.5  mm2).  Con¬ 
sidering  an  interrogation  overlap  of  47  percent  means 
that  about  18  vectors  would  be  calculated  within  the 
vortex  core.  Given  that  the  size  of  the  vortex  grows 
as  it  ages,  the  number  of  vectors  calculated  within 
the  core  would  also  increase.  In  this  present  test,  for 
example,  a  vortex  that  had  grown  to  a  diameter  of 
20  mm  would  contain  about  72  vectors. 

A  single  file  for  each  wake  age  included  a  flag 
for  every  velocity  measurement  that  indicated  the 
nature  of  that  vector  (such  as  valid,  invalid,  interpo¬ 
lated,  recalculated,  or  not  calculated).  Reduced  data 
files  containing  the  coordinates  of  each  calculation 
node  and  the  three  components  of  velocity  (in  terms 
of  displacement)  were  stored  in  ASCII  format  and 
then  transferred  to  a  mainframe  computer  for  anal¬ 
ysis.  The  velocity  components  were  first  converted 
from  displacement  units  to  velocity  units  based  on 
the  pulse  duration  for  that  particular  measurement. 
The  physical  coordinates  were  transformed  so  that 
r  =  z  =  0  would  correspond  to  the  center  of  the 
rotor  hub  (fig.  2). 

Because  of  centripetal  effects,  heavier  seed  par¬ 
ticles  appear  to  have  been  driven  outward  from  the 
center  of  the  vortex  as  it  moved  through  the  fluid. 
Because  the  size  of  the  remaining  lighter  particles 
were  most  likely  much  smaller,  less  light  would  have 
been  scattered  from  this  region,  thereby  giving  it  the 
appearance  of  a  “void.”  Nevertheless,  the  particles 
that  were  present  inside  the  void  were  sufficient  to 
obtain  PIV  measurements.  In  searching  for  vortex 
characteristics  that  can  be  known  with  certainty,  the 
size  of  the  “apparent  void”  cannot  always  be  equated 
to  the  size  of  the  core,  nor  can  the  center  of  the 
void  be  assumed  to  be  coincident  with  the  center  of 
swirl.  On  the  other  hand,  the  physical  location  of 


the  void  was  an  indisputable  feature  and  was,  there¬ 
fore,  judged  to  be  the  most  reliable  indicator  of  how 
the  vortex  was  convecting  through  the  fluid.  In  addi¬ 
tion,  in  order  to  account  for  the  ever-present  vortex 
meander,  the  coordinates  of  the  voids  for  a  set  of 
images  acquired  over  time  at  a  particular  azimuth 
behind  the  rotor  provided  a  very  accurate  basis  for 
conditionally  averaging  the  images. 

The  distribution  of  particles  around  an  apparent 
vortex  void  often  exhibited  a  peculiar  pattern.  Parti¬ 
cles  were  not  always  uniformly  distributed  in  the  flow 
as  it  approached  the  rotor.  In  these  instances,  some 
form  of  particle  “banding”  frequently  occurred  adja¬ 
cent  to  the  void  region.  Calculations  have  been  made 
that  show  that  banding  is  also  a  natural  consequence 
of  the  balance  of  forces  (drag,  centrifugal,  and  Cori¬ 
olis)  that  act  on  the  seed  particles  within  a  vortex 
(ref.  20).  Small  particles  settle  in  an  annular  region 
that  is  smaller  than  the  vortex  core,  compared  to 
larger  particles  that  yield  a  band  that  is  larger  than 
the  vortex  core.  Regardless  of  how  uniform  the  flow 
appeared  to  be  seeded  in  the  present  experiment,  a 
fairly  broad  range  of  void  sizes  was  observed  at  all 
wake  ages  (fig.  3).  The  average  size  of  the  void  in¬ 
creased  rapidly  during  the  first  100°  of  azimuth,  after 
which  the  rate  of  increase  diminished  substantially. 

In  this  study,  the  sequence  for  extracting  the 
characteristics  of  the  trailing  vortex  began  with  a 
physical  examination  of  each  raw  image  to  deter¬ 
mine  the  size  and  location  (relative  to  a  superim¬ 
posed  mesh)  of  the  apparent  void  (fig.  4).  At  each 
wake  age  the  images  were  averaged  to  smooth  out 
the  small  irregularities  that  were  invariably  present 
in  the  flow.  Because  the  vortex  was  located  at  dif¬ 
ferent  positions  from  image  to  image  (vortex  mean¬ 
der)  ,  performing  a  simple  average  would  have  risked 
smearing  out  important  details  concerning  the  struc¬ 
ture  of  the  vortex.  Instead,  the  data  were  averaged 
only  after  aligning  the  set  of  images  based  on  the 
most  dominant  recurring  feature  in  the  flow,  which 
in  this  case  was  the  apparent  void. 

This  procedure  is  referred  to  as  conditional  en¬ 
semble  averaging,  the  condition  here  being  the  align¬ 
ment  of  all  the  images  in  the  set  based  on  their  void 
centers.  The  image  with  a  void  center  closest  to  the 
mean  location  that  was  established  for  the  set  of  im¬ 
ages  was  selected  as  the  “anchor.”  The  indices  of 
all  image  matrices  (that  is,  grid  points  i,j)  for  each 
wake  age  were  then  adjusted  according  to  the  offset 
of  each  void  from  the  coordinates  of  the  void  in  the 
anchor  image. 

At  this  intermediate  point,  the  conditionally  av¬ 
eraged  flow  field  at  each  wake  age  embodies  the  trail¬ 
ing  vortex  within  a  highly  nonuniform  downwash. 
The  time-averaged  downwash  for  a  single-bladed  ro¬ 


tor  can  be  determined  by  averaging  the  flow  field 
over  an  entire  revolution  of  the  rotor  (refs.  5,  7). 
When  the  time- averaged  downwash  velocity  field  is 
subtracted  from  the  phase-averaged  velocity  field, 
the  resulting  velocity  profiles  tend  to  appear  more 
symmetric  along  the  extremities  of  the  vortex,  but 
not  near  the  core  radius  (ref.  5).  In  the  present  ex¬ 
periment,  this  procedure  could  not  be  implemented 
because  data  were  not  obtained  over  the  required 
half  revolution  in  each  zone.  However,  the  axial  and 
radial  locations  of  the  voids  at  each  wake  age,  along 
with  the  known  time  intervals  between  each  wake 
age,  were  known  and  did  yield  the  convection  ve¬ 
locity  of  the  vortex  at  each  wake  age.  These  com¬ 
ponents  of  the  convection  velocity,  uc  and  wc,  were 
subtracted  from  the  phase-averaged  flow  field.  How¬ 
ever,  because  the  axial  component  of  the  downwash 
over  the  inboard  region  was  significantly  higher  than 
the  axial  convection  velocity  of  the  vortex  (u>c),  ve¬ 
locity  profiles  across  the  center  of  the  vortex  cannot 
be  expected  to  be  symmetric. 

To  extract  information  about  the  geometry  of 
the  vortex  (after  performing  averages  based  on  the 
void  centers  at  each  wake  age),  the  center  of  the 
vortex  was  assumed  to  be  at  the  center  of  swirl2. 
Using  the  void  center  as  a  starting  point,  the  sur¬ 
rounding  locations  were  interrogated  to  determine 
the  best  node  for  which  the  sum  of  the  dot  prod¬ 
ucts  of  two  unit  vectors  produced  a  minimum  over 
a  neighborhood3  of  locations  surrounding  the  candi¬ 
date  node.  The  coordinates  of  the  neighboring  node 
relative  to  the  candidate  node  defined  one  of  the  unit 
vectors.  The  velocity  vector  at  the  neighboring  node 
defined  the  other  unit  vector.  With  the  angles  of 
these  two  unit  vectors  denoted  by  a  and  /?,  the  node 
nearest  the  center  of  swirl  was  said  to  be  found  when 
the  following  expression  was  a  minimum: 

N 

]T|cos(o  — /?)|/A  (1) 

n=l 

where  N  is  the  total  number  of  neighboring  points 
considered.  Relative  to  the  center  of  swirl  for  a  pure 
vortex  in  a  stationary  flow,  a  and  0  will  be  orthog¬ 
onal  for  all  neighboring  points  and  the  sum  will  be 
exactly  zero.  For  most  real  flows  the  sum  will  not  be 
precisely  zero,  and  the  procedure  works  best  when 
the  neighborhood  under  consideration  does  not  ex¬ 
tend  beyond  about  one  or  two  core  diameters  (which 
is  normally  sufficient  to  cover  most  of  the  vorticity 

2  Neither  the  void  center  nor  the  swirl  center  was  necessarily 
at  the  center  of  the  vorticity  that  was  accreted  from  the  rotor 
wake,  itself  having  the  form  of  a  spiraling  sheet  containing 
varying  amounts  of  vorticity. 

3Because  the  quality  of  the  data  within  the  void  may  be  in 
question,  these  vectors  were  omitted  during  this  operation. 


in  a  single  trailing  vortex) . 

The  size  of  the  vortex  core  was  found  by  first  di¬ 
viding  the  surrounding  region  into  annular  zones  and 
calculating  the  average  swirl  velocity  and  the  associ¬ 
ated  mean  radius  of  each  annulus.  These  data  were 
fit  with  a  least-squares  spline  subject  to  the  condi¬ 
tion  that  the  resulting  curve  have  only  one  inflection 
along  its  inner  extent  (nominally  set  to  twice  the  esti¬ 
mated  core  radius) .  The  point  of  inflection  was  used 
to  define  the  radius  of  the  vortex  core,  because  that 
is  where  a  maximum  value  for  the  swirl  velocity  is 
reached.  This  procedure  was  considered  to  offer  a 
more  rational  approach  for  determining  the  size  of 
the  vortex,  especially  when  the  velocity  peaks  that 
are  characteristic  of  vortices  (and  upon  which  core 
sizes  are  traditionally  based)  are  dependent  on  how 
the  vortex  is  sliced  (vertical,  horizontal,  or  other¬ 
wise)  . 

Vorticity  was  calculated  by  integration  (circula¬ 
tion  box  method)  according  to: 

uie  =  A~1T  =  A~l  j>  (u,w)  ■  d{f,z)  (2) 

where  the  direction  of  integration  was  such  that  the 
enclosed  area,  A,  was  on  the  left  of  the  integration 
path.  The  size  of  each  circulation  box  was  defined 
by  the  local  grid  spacing,  which  on  average  rendered 
A  =  2.5  mm2. 

TEST  CONDITIONS 

The  experiment  was  performed  at  a  constant  ro¬ 
tor  speed  of  870  rpm  (14.5  Hz),  which  for  a  rotor  di¬ 
ameter  of  191.5  cm  corresponds  to  a  tip  speed  of  87.2 
m/s.  Based  on  an  average  ambient  temperature  of 
16.7  degrees  Centigrade  (°C),  the  Reynolds  number 
based  on  rotor-tip  chord  was  Re  =  0.6  x  106  and  the 
Mach  number  was  M  =  0.26.  The  collective  pitch 
angle  of  the  rotor  was  fixed  at  11°. 

To  obtain  high- resolution  PIV  images,  the  cam¬ 
eras  were  focused  on  a  small  region  of  the  wake  be¬ 
hind  the  trailing  edge  of  the  rotor.  The  camera  view 
was  radially  and  axially  translated  to  provide  three 
zones  of  coverage  (A,  E,  and  F).  The  laser  pulse  de¬ 
lays  needed  to  capture  the  highest  velocities  ranged 
from  3.5  to  5.0  fis.  Results  presented  later  in  this  re¬ 
port  show  that  even  for  the  longest  pulse  delay,  5  mi¬ 
croseconds  (/us) ,  the  vortex  moved  less  than  1/1000  of 
the  core  radius  between  pulses.  The  blade  was  con¬ 
sidered  to  be  at  zero  azimuth  when  the  trailing  edge 
of  the  rotor  blade  was  horizontal.  When  the  blade  is 
in  this  position,  the  light  sheets  (both  of  which  in¬ 
tersect  in  a  horizontal  plane)  pass  just  beneath  the 
trailing  edge  of  the  blade.  The  procedure  for  taking 
images  at  the  different  wake  ages  was  accomplished 


by  delaying  the  1/rev  signal  from  the  rotor,  thereby 
allowing  the  tip  of  the  rotor  blade  to  move  past  the 
measurement  plane  by  the  desired  azimuthal  incre¬ 
ment.  The  data  were  taken  at  a  nominal  increment 
of  15°  (fig.  5),  with  some  exceptions  as  required  to 
avoid  reflections  of  the  laser  light  from  any  portion 
of  the  blade  surface  that  might  damage  the  camera 
sensors. 

Rotor  loads  and  tunnel  velocity  measurements 
(table  1)  were  taken  from  an  earlier  experiment  (ref. 
15).  Strictly  speaking,  the  rotor  is  not  in  hover,  as 
evidenced  by  the  magnitude  of  the  velocity  upstream 
of  the  rotor.  Because  the  flow  around  the  rotor  is 
bounded  by  the  wind  tunnel  walls,  and  because  the 
rotor  wake  lowers  the  pressure  in  the  test  section  (in 
turn  inducing  a  finite  flow  everywhere  ahead  of  the 
rotor),  a  slight  climb  condition  cannot  be  avoided 
(see  ref.  15  for  a  full  discussion).  Assuming  the  up¬ 
stream  velocity  entirely  represents  the  climb  velocity, 
then  the  climb  rate  for  this  part  of  the  experiment 
can  be  calculated  from  wrj0/ R^l. 


Table  1:  Rotor  loads 


Rotational  speed  of  rotor,  rpm 

870 

Collective  pitch  angle,  deg 

11 

Solidity 

0.0691 

Torque,  foot-pounds  (ft-lbs) 

11.69 

Thrust,  lb 

39.76 

Upstream  velocity,  ft/s 

1.28 

Climb  rate 

0.0045 

Hover  figure  of  merit 

0.61 

CT 

0.00657 

Cq 

0.00062 

Ct  /  cr 

0.095 

RESULTS  AND  DISCUSSION 

The  center  of  the  vortex  void,  the  center  of  vor¬ 
ticity,  the  center  of  swirl,  and  the  center  between 
velocity  peaks  are  all  considered  to  be  distinct  quan¬ 
tities.  For  each  wake  age,  starting  with  =  0°,  the 
flow  field  has  been  summarized  in  terms  of  quantities 
that  focus  on  the  trailing  vortex  (he.,  fig.  6).  Vortex 
wander  was  qualitatively  identified  by  plotting  the 
locations  of  the  vortex  voids  for  each  image  set  (up¬ 
per  left  corner  of  figure) .  To  the  right  of  this  plot  is  a 
sketch  of  the  rotor  showing  the  location  of  the  rectan¬ 
gular  zone  covered  by  the  entire  calculation  grid  (92 
steps  in  the  z  direction  and  79  steps  in  the  r  direc¬ 
tion).  The  three  field  plots  in  the  left- most  column 
show  the  out-of-plane  vorticity,  (contours  with  red 
and  blue  hues  to  indicate  direction  and  a  continuous 
range  of  saturations  to  indicate  magnitude),  inplane 


velocity  relative  to  the  convection  velocity  of  the  vor¬ 
tex,  ur  +  wr  (vectors  with  a  continuous  range  of  hues 
to  indicate  magnitude),  and  the  out-of-plane  veloc¬ 
ity,  v  (contours  with  red  and  blue  hues  to  indicate 
direction  and  a  continuous  range  of  saturations  to 
indicate  magnitude).  The  three  profiles  in  the  cen¬ 
ter  column  show  ug,  wr,  and  v  along  z  —  Constant 
and  passing  through  the  center  of  swirl.  The  central 
portion  of  these  profiles  is  shown  in  red  and  demar¬ 
cated  with  dashed  vertical  lines  in  order  to  identify 
the  void  region.  The  horizontal  and  vertical  dashed 
lines  included  in  each  of  the  field  and  profile  plots 
mark  the  location  of  the  swirl  center. 

The  lower  plot  in  the  right-most  column  shows 
the  inplane  swirl  velocity,  ip,  averaged  over  incremen¬ 
tal  annular  regions  with  radii  (£)  measured  from  the 
swirl  center.  The  upper  plot  in  the  right-most  col¬ 
umn  shows  an  enlargement  of  the  inplane  velocity 
(ur  +  wr),  limited  to  a  region  twice  the  diameter  of 
the  vortex  core.  The  two  circles  identify  the  size 
of  the  vortex  core  (solid  curve)  and  the  vortex  void 
(dashed  curve).  Finally,  values  are  given  at  each 
wake  age  for  the  void  center,  the  vorticity  center 
(or  centroid  of  vorticity),  the  swirl  center,  the  vor¬ 
tex  convection  velocity  (obtained  from  the  path  of 
the  vortex  voids) ,  and  the  core  radius  along  with  the 
maximum  swirl  velocity  (both  of  which  are  fixed  by 
the  peak  value  found  from  the  curve  fit  to  <p  vs.  £). 
These  values  are  also  given  in  tables  2  to  4. 

Based  on  a  review  of  the  data  for  all  wake  ages 
(selected  cases  are  given  in  figs.  6-10),  the  highest 
velocities  are  found  to  occur  within  the  first  15  de¬ 
grees  behind  the  rotor  blade.  During  this  interval, 
the  maximum  inplane  (swirl)  and  out-of-plane  ve¬ 
locities  reached  44  and  12  percent  of  the  rotor  tip 
velocity  (87.2  m/s),  respectively.  Although  the  de¬ 
gree  of  vortex  wander  was  somewhat  erratic,  it  gen¬ 
erally  increased  with  wake  age.  On  average,  the 
standard  deviation  (based  on  the  mean  location  of 
the  void  center)  was  only  about  9  percent  of  the  di¬ 
ameter  of  the  vortex  core  (the  worst  case  being  18 
percent).  The  calculated  locations  of  the  vorticity, 
swirl,  and  velocity-peak4  centers  generally  differed 
from  the  void  center  by  about  11,  26,  and  11  percent 
of  the  vortex  diameter,  respectively  (fig.  12).  There 
was  no  apparent  correlation  between  the  void  and 
swirl  centers  from  image  to  image. 

The  out-of-plane  velocity  exhibited  the  most  un¬ 
usual  appearance.  Aside  from  the  boundary  layer 
wake  that  is  shed  from  the  trailing  edge  of  the  ro¬ 
tor  blade  and  is  most  distinct  when  ip  <  30° ,  there 
appear  to  be  two  concentrated  zones  of  oppositely 
directed  flow  near  the  center  of  the  vortex.  Because 

4  The  velocity-peak  center  is  based  on  the  locations  that 
define  df,  (see  fig.  11). 


both  of  these  zones  are  generally  offset  from  the  vor¬ 
tex  center,  the  bisecting  profiles  taken  along  z  = 
Constant  do  not  fully  reveal  the  magnitude  of  the 
“excess”  and  “deficit”  velocities  that  contribute  to 
the  overall  structure  of  the  trailing  vortex. 

Even  though  the  convection  velocity  of  the  vor¬ 
tex  at  each  wake  age  has  been  subtracted,  the  rela¬ 
tive  “vertical”  component  of  velocity,  wr ,  exhibits  an 
asymmetric  profile.  Greater  symmetry  near  the  vor¬ 
tex  may  be  achieved  by  subtracting  the  local  down- 
wash,  which  can  only  be  approximated  from  a  time 
average  of  the  flow  field  at  corresponding  locations 
over  one  or  more  revolutions  (refs.  5,  7).  However, 
it  is  not  clear  that  this  effort  will  yield  any  better 
understanding  of  the  vortex  structure,  especially  be¬ 
cause  the  resulting  profiles  (including  the  peaks)  will 
still  be  asymmetric  (refs.  5,10).  Given  the  inevitable 
likelihood  of  asymmetry,  traditional  estimates  of  the 
core  size  and  peak  velocity  would  still  have  to  be 
based  on  an  average  between  the  two  halves  of  the 
asymmetric  profile,  and  this  is  effectively  what  is  al¬ 
ready  provided  in  the  more  general  annular  averages 
that  produce  the  <p  vs.  £  plot.  Nevertheless,  the  pro¬ 
file  plots  are  important  because  they  retain  certain 
nuances  that  can  prove  helpful  in  explaining  other 
unexpected  features  in  the  flow. 

During  the  initial  development  of  the  trailing 
vortex,  ip  <  30°,  the  vorticity  within  the  vortex 
reached  a  peak  value  very  near  the  center  of  the  vor¬ 
tex,  and  rapidly  diminished  to  near  zero  at  a  distance 
not  far  beyond  the  core  diameter.  For  larger  wake 
ages,  ip  >  45°,  the  peak  values  of  vorticity  appear  to 
reside  in  an  annular  region  defined  by  the  size  of  the 
void,  yet  are  clearly  within  the  core  diameter.  As¬ 
suming  these  measurements  are  valid,  this  variation 
in  vorticity  would  imply  that  much  of  the  interior 
region  of  the  vortex  was  not  fully  in  solid-body  rota¬ 
tion.  Evidence  for  this  can  be  found  in  the  vertical 
component  of  velocity,  wr ,  which  shows  a  change  in 
the  slope  near  the  vortex  center  (see,  for  example, 
the  profile  at  ip  =  195°  appearing  in  fig.  9).  This 
peculiar  departure  from  the  traditional  view  of  a  vor¬ 
tex  has  appeared  in  other  rotor  vortex  data  (ref.  6, 
where  the  flow  was  seeded  with  incense  particles  and 
measurements  made  with  a  laser  Doppler  velocime- 
try  (LDV)  technique,  and  ref.  24). 

To  be  sure,  other  explanations  (seeding  being 
among  them)  for  the  unusual  velocity  and  vorticity 
profiles  over  the  interior  of  the  vortex  are  plausible. 
Although  centripetal  effects  may  have  been  a  factor, 
the  extent  to  which  the  velocity  within  the  core  was 
significantly  altered  is  unknown  because  the  size  of 
the  particles  was  not  measured.  Another  factor  may 
have  been  the  rival  between  legitimate,  but  barely 
visible,  small  particles  and  the  presence  of  sensor 


noise,  the  result  being  an  abundance  of  interpolated 
vectors  within  the  void  region.  Until  these  issues 
can  be  more  thoroughly  addressed,  it  is  uncertain 
whether  the  data  show  an  onset  of  vortex  instability; 
reflect  a  difference  in  particle  dynamics  because  of  a 
nonuniform  distribution  of  particle  sizes  in  a  highly 
accelerated  flow;  or  are  simply  due  to  a  processing 
limitation  arising  from  a  marginal  scattering  of  light 
from  very  small  particles  that  resulted  in  incorrect 
interpolations. 

The  location  with  the  greatest  number  of  in¬ 
terpolations  always  occurred  within  the  void  region. 
The  number  of  interpolated  values  contributing  to 
the  average  at  these  “worst-case”  locations  ranged 
from  16  to  60  percent.  Although  far  from  ideal,  no 
location  consisted  solely  of  interpolated  values. 

To  address  a  possible  flaw  in  the  logic  used  to 
average  the  data  during  the  analysis  phase,  the  case 
for  ip  =  195°  was  reexamined.  First,  the  results  for 
single  images  (hence,  unaveraged)  were  individually 
studied  to  see  if  the  nonlinear  slope  in  the  core  ve¬ 
locity  might  simply  have  been  an  artifact  of  the  av¬ 
eraging  process.  Although  the  unaveraged  velocity 
profiles  were  noticeably  more  “noisy,”  the  peculiar 
change  in  slope  near  the  center  of  the  vortex  was 
found  to  be  present  and  similar  to  the  conditionally 
averaged  results.  As  a  final  check,  a  different  cri¬ 
terion  was  considered  for  obtaining  the  conditional 
average.  In  this  case  the  more  classical  approach  of 
using  the  midpoint  between  the  velocity  peaks  (high¬ 
est  maximum  and  minimum  values  occurring  dur¬ 
ing  horizontal  cuts  across  the  vortex)  was  adopted. 
Again,  the  variation  in  slope  across  the  vortex  was 
essentially  unchanged  from  that  observed  when  the 
conditional  averaging  criterion  was  based  on  the  lo¬ 
cations  of  the  voids.  In  fact,  the  same  conclusion  was 
reached  after  considering  the  swirl  centers  as  the  av¬ 
eraging  criterion. 

Another  curious  feature  of  the  vortex  is  its 
slightly  elliptical  shape.  This  was  also  observed  in 
the  flow  pattern  constructed  from  LDV  measure¬ 
ments  in  the  near  wake  of  a  hovering  rotor  (ref.  22)  A 
minor  amount  of  ellipticity  can  be  attributed  to  the 
fact  that  the  axis  of  the  vortex  is  not  perpendicular 
to  the  plane  of  the  light  sheet.  This  would  cause  the 
vortex  to  appear  stretched  in  the  radial  and  axial 
directions  (relative  to  the  rotor  coordinates).  The 
amount  of  stretching  in  both  directions  was  nearly 
the  same  during  the  early  wake  ages,  thereby  in¬ 
creasing  the  apparent  core  size  by  about  0.06  per¬ 
cent  and  having  essentially  no  effect  on  ellipticity. 
As  the  vortex  aged,  the  apparent  stretching  became 
progressively  more  biased  in  the  axial  direction  (be¬ 
cause  the  radial  displacement  was  approaching  a  con¬ 
stant)  ,  but  the  degree  of  ellipticity  caused  by  the  an¬ 


gle  of  the  vortex  axis  was  still  less  than  0.4  percent 
after  one  revolution.  Judging  by  the  amount  of  ellip¬ 
ticity  observed  in  the  data  (sometimes  reaching  an 
aspect  ratio  of  about  1.23),  and  by  the  alternating 
orientation  of  the  major  axis,  the  cause  for  this  ellip¬ 
tical  flow  pattern  cannot  be  due  to  a  changing  angle 
of  the  vortex  axis  as  the  vortex  ages  (unless  small- 
amplitude  undulations  are  superimposed  on  the  vor¬ 
tex  trajectory,  although  the  data  in  figure  12  do  not 
show  this). 

To  study  this  phenomenon  in  more  detail,  the 
dimensions  of  the  vortex  (major  and  minor  axes) 
were  equated  to  the  distances  between  the  velocity 
peaks  found  from  horizontal  and  vertical  cuts  across 
the  center  of  the  vortex  (recall  fig.  11).  The  peak- 
to-peak  distances  associated  with  the  horizontal  and 
vertical  cuts  were  designated  by  dh  and  dv ,  respec¬ 
tively,  and  these  are  shown  for  wake  ages  ip  =  0°  to 
390°.  Using  red  and  green  color  bars  to  indicate  the 
size  of  dh  relative  to  dV)  the  data  suggest  that  the 
elliptical  pattern  is  rotating.  This  figure  also  shows 
that  the  nonlinear  behavior  of  the  velocity  profile 
within  the  core  does  not  depend  on  the  direction  of 
the  cut  (horizontal  or  vertical)  or  the  direction  of  the 
major  axis  of  the  ellipse  (horizontal  or  vertical). 

According  to  the  location  of  the  void  center,  the 
radial  and  axial  coordinates  of  the  vortex  with  in¬ 
creasing  wake  age  are  known  (fig.  13).  Based  on 
curve  fits  to  this  data,  the  corresponding  convection 
velocities  were  calculated.  These  data  indicate  that 
the  highest  radial  velocity  occurred  immediately  af¬ 
ter  the  trailing  vortex  left  the  trailing  edge  of  the  ro¬ 
tor.  This  figure  shows  that  the  magnitude  of  the  ra¬ 
dial  velocity  rapidly  decreased  and  then  approached 
an  asymptotically  diminishing  value  near  the  com¬ 
pletion  of  one  rotor  revolution.  The  behavior  of  the 
axial  velocity  component  was  quite  different.  The 
magnitude  of  the  axial  velocity  was  lowest  imme¬ 
diately  after  the  trailing  vortex  left  the  rotor,  but 
then  increased  until  its  progress  was  noticeably  im¬ 
peded  after  the  passage  of  the  following  rotor  blade. 
The  velocities  may  have  a  slightly  different  appear¬ 
ance  if  curve  fits  to  the  data  are  separated  into  seg¬ 
ments  before  and  after  passage  of  the  second  blade 
(ip  =  180°).  The  wake  contraction  can  also  be  viewed 
in  the  r-z  plane  (fig.  12),  which  shows  the  character¬ 
istic  “knee”  due  to  the  influence  of  the  second  blade, 
and  the  increasing  dominance  of  the  axial  velocity  as 
the  wake  aged. 

Based  on  the  radial  distribution  of  the  swirl  ve¬ 
locity  at  each  wake  age,  the  peak  (or  inflection)  in 
the  curve  fit  was  used  to  define  the  core  size  and 
the  maximum  swirl  velocity.  The  initial  diameter  of 
the  vortex  core  was  about.  10  percent  of  the  chord  of 
the  rotor  blade  (c  =  104  mm)  and  more  than  dou- 


bled  its  size  after  one  revolution  of  the  rotor  (fig. 
12).  The  maximum  swirl  velocity  was  highest  during 
the  initial  formation  of  the  trailing  vortex,  decreased 
rapidly  during  the  first  half  revolution,  and  settled 
on  a  nearly  constant  value  (on  average)  during  the 
last  half  revolution  of  the  rotor.  It  is  difficult  to  say 
exactly  how  the  intensity  of  the  vortex  (which  is  re¬ 
lated  to  the  maximum  swirl  velocity)  changed  with 
wake  age  because  of  the  unavoidable  influence  of  the 
downwash  on  the  magnitude  of  the  swirl  velocity  (the 
vortex  is  embedded  in  a  nonuniform  velocity  field) . 

To  examine  more  closely  the  structure  of  the 
trailing  vortex,  contour  plots  of  the  vertical  compo¬ 
nent  of  velocity,  wr,  and  the  out-of-plane  vorticity, 
tig ,  are  shown  at  wake  ages  ip  =  15°  and  195°  along 
with  profiles  taken  from  several  closely  spaced  hori¬ 
zontal  cuts  across  the  vortex  (figs.  14).  The  velocity 
profiles  are  clearly  more  linear  within  the  core  of  the 
vortex  at  ip  =  15°  than  they  are  at  ip  =  195°.  The 
corresponding  vorticity  profiles  exhibit  a  single  peak 
when  crossing  the  vortex  at  ip  =  15° ,  whereas  a  dou¬ 
ble  peak  is  evident  along  several  cross  sections  of  the 
flow  at  ip  —  195°. 

Another  point  of  interest  is  the  radial  distribu¬ 
tion  of  circulation  across  the  vortex.  Using  the  swirl 
velocity  to  calculate  the  circulation,  and  normalizing 
both  <p  and  F  by  values  existing  at  the  core  radius 
(CO,  results  for  ip  —  15°,  195°,  and  375°  are  shown  in 
figure  15.  As  the  distance  increased  from  the  center 
of  swirl  (C  =  0) ,  the  circulation  rapidly  increased  un¬ 
til  reaching  a  point  just  beyond  the  core  radius,  after 
which  the  rate  of  growth  decreased  significantly. 

Similar  to  the  selected  views  of  circulation,  the 
radial  distribution  of  vorticity  was  examined  during 
the  development  of  the  trailing  vortex.  However,  be¬ 
cause  the  vorticity  center  and  the  swirl  center  (as 
well  as  the  void  center)  data  were  generally  not  collo¬ 
cated,  the  normalized  vorticity  profile  was  referenced 
to  the  vorticity  center.  The  peak  vorticity  occurred 
at  the  center  ( k  =  0)  when  ip  =  15°  (fig.  16).  At 
subsequent  wake  ages  (ip  =  195°  and  375°),  the  vor¬ 
ticity  in  the  center  decreased  and  the  peaks  progres¬ 
sively  shifted  toward  the  core  diameter.  A  similar 
change  in  the  distribution  of  vorticity  was  reported 
in  an  earlier  fixed-wing  experimental  study  using  a 
vorticity  meter  (ref.  23).  As  discussed  earlier,  it  is 
not  known  to  what  extent  this  characteristic  might 
be  due  to  potentially  inaccurate  data  within  the  void 
region.  Nevertheless,  the  fact  that  this  radial  shift 
in  the  vorticity  peak  has  been  observed  while  using 
other  measurement  techniques  (vorticity  meter  and 
LDV)  suggests  that  the  PIV  data  are  authentic  and 
that  the  vortex  was  indeed  becoming  unstable. 

The  self-similar  nature  of  the  swirl  velocity  and 
the  vorticity  for  the  entire  set  of  wake  ages  (ip  =  0° 


to  390°)  are  shown  in  figure  17.  The  swirl  veloc¬ 
ity  was  made  dimensionless  by  its  value  at  the  core 
radius  (<p/ipc),  and  the  radius  (measured  from  the 
swirl  center)  was  made  dimensionless  by  the  core  ra¬ 
dius  (C/Cc)-  Some  data  points  for  the  dimensionless 
swirl,  Tp,  exceed  1.0  because  the  maximum  swirl  ve¬ 
locity,  Pc,  was  based  on  curve  fits  to  <p  vs.  tp  at  each 
wake  age.  The  dimensionless  form  of  a  popular  vor¬ 
tex  model  (ref.  18) 


has  also  been  included  in  the  plot.  The  best  fit  to  the 
data  occurred  when  n  =  3  and  r*/rc  =  1.26.  Ac¬ 
cording  to  this  model,  the  circulation  around  a  con¬ 
tour  placed  at  the  core  boundary,  rc,  was  about  79 
percent  of  the  total  far-field  (or  large-radius)  value, 
T*.  Note  that  at  the  core  radius  Tp  —  tp  =  1,  so  that 
equation  (3)  requires  that  r*/rc  =  21/"  =  1.26  when 
n  —  3. 

The  self-similar  distribution  of  accumulated  vor¬ 
ticity  can  be  found  from  the  product  of  the  local 
vorticity  (ojg)  and  the  incremental  area  over  which 
it  acts  (ArAz),  and  then  summed  over  regions  de¬ 
fined  by  the  radius  from  the  vorticity  center  ( k ).  In 
principle,  the  result  should  be  equivalent  to  the  vari¬ 
ation  of  circulation  with  radius,  provided  the  swirl 
and  vorticity  centers  are  coincident.  Within  the  core 
( k  <  pc),  the  distribution  of  vorticity  is  quite  similar 
and  is  well  represented  by  the  model,  tpp.  Beyond 
the  core  radius  the  accumulated  vorticity  asymp¬ 
totically  approaches  a  constant  value  and  on  aver¬ 
age  agrees  with  that  predicted  by  the  vortex  model, 
which  states  that  Tpp  — >■  1.26  as  p  becomes  large. 

The  core-radius  circulation,  rc,  was  observed  to 
increase  slightly  with  wake  age  (fig.  18),  Because 
some  vorticity  was  found  to  exist  beyond  the  core 
radius  at  all  wake  ages,  and  because  the  core  ra¬ 
dius  generally  increased  with  wake  age  (fig.  18),  it 
is  probable  that  additional  vorticity  was  gradually 
being  gathered  into  the  core.  The  total  far-field  vor¬ 
ticity  appears  to  be  fairly  constant,  F»  =  1.42  m2/s. 
The  growth  in  the  size  of  the  vortex  core  can  be 
approximated  by  the  following  model  (ref.  24,  and 
further  developed  in  ref.  22): 

pc  =  (Pt0  +  Mvip/n)1'2  (4) 

where  S  =  16  provides  the  best  fit  to  the  present 
data  (for  laminar  flow,  S  =  1).  Data  included  from 
another  experiment  (ref.  22,  based  on  LDV  measure¬ 
ments)  was  reported  to  be  best  fit  by  the  model  when 
J  =  8.  It  should  be  noted  that  <5  strongly  depends 
on  the  initial  core  size,  CcO,  and  that  lower  values  of 


(cO  require  lower  values  of  S.  The  results  presented 
in  reference  25,  which  are  for  a  single-bladed  rotor, 
show  a  similar  range  of  core  sizes  and  growth  rate. 

Based  on  an  isolated  vortex  model  (refs.  20, 
25)  that  assumes  a  turbulence  mixing  length  that 
is  proportional  to  the  vortex  radius,  the  ratio  of  the 
core-radius  to  large-radius  circulation  should  be  fixed 
by  the  vortex  Reynolds  number,  T*/v.  However,  the 
data  for  this  test  (T»/u  ss  1.2  x  105)  fall  above  the 
curve  and  are,  on  average,  aligned  with  a  laminar 
value  for  rc/r*  (fig.  18).  This  would  imply  that 
S  ss  1  rather  than  16  as  required  for  a  good  fit  to 
the  (c  vs.  data,  and  would  suggest  that  either  the 
models  are  imprecise  or  that  the  core  radius  is  not 
as  large  as  measured. 

Just  as  the  core-radius  circulation  was  less  than 
the  large-radius  circulation  during  the  first  revolu¬ 
tion  of  the  rotor,  the  large-radius  circulation  is  ex¬ 
pected  to  be  less  than  the  maximum  bound  circu¬ 
lation  on  the  blade.  Although  the  blade  circulation 
was  not  measured  in  this  test,  a  value  can  be  ap¬ 
proximated  from  blade-element  momentum  theory 
(ref.  26).  For  a  blade  with  constant  chord,  the  to¬ 
tal  thrust  coefficient  can  be  stated  in  terms  of  the 
dimensionless  inflow,  A,  as 

Ct  =  f  4A 2rdr  (5) 

Jo 


expressions  for  the  inflow  in  the  form  of  A  =  Aor" 
have  been  proposed  (refs.  20,  26).  The  results  for 
n  =  0  (constant  inflow),  n  =  0.5  (linear  circulation), 
and  n  =  1  (linear  inflow)  are  shown  in  figure  18. 
Although  this  power  function  has  served  well  to  il¬ 
lustrate  various  performance  trends,  the  inflow  used 
in  the  present  analysis  was 


A  =  A0y  sin  ^-r  (9) 

From  equations  (5),  (7),  and  (9)  the  maximum  bound 
circulation  on  the  rotor  blade  is 


(10) 


Recalling  that  Ct  =  0.00657  and  a  =  0.0691,  then 
equation  (10)  states  that  (F{,)ma:E  =  2.13  m2/s  for 
this  test.  To  summarize: 


^  =  0.79  (11) 

-A- =  0.67  (12) 

A  bmax 

implying  that  the  core-radius  circulation  is  only 
about  53  percent  of  the  maximum  bound  circulation 
on  the  rotor  blade. 


and  the  sectional  value  in  terms  of  the  local  circula¬ 
tion,  Ti,,  on  the  blade  as 


dC-  =  <6> 

From  the  above  two  equations  it  follows  that  the 
inflow  velocity  and  the  blade  circulation  (both  are 
functions  of  r)  are  related  by 


QRc 


=  4A2 


(7) 


The  local  inflow  velocity  for  hover  can  be  stated  as 
(ref.  26) 


_  crCig 
16 


1  + 


32  rO 

&Cla 


(8) 


where  9  is  the  local  blade  angle.  Because  the  climb 
rate  in  this  test  was  insignificant,  hover  conditions 
are  assumed  to  apply.  The  change  in  blade  angle 
along  the  span  for  the  rotor  used  in  this  test  is  shown 
in  figure  18.  If  the  blade  angle,  9(r),  is  known,  and  if 
it  is  assumed  that  Cia  &  0.1,  the  local  inflow  veloc¬ 
ity  predicted  by  equation  (8)  can  be  determined  (fig. 
18).  However,  in  order  to  more  easily  perform  the 
integration  given  in  equation  (5),  simple  analytical 


CONCLUSIONS 

1 .  The  highest  velocities  occurred  within  the  first 
15°  behind  the  rotor  blade.  The  swirl  and  axial  ve¬ 
locities  reached  44  and  12  percent  of  the  rotor  tip 
speed,  respectively. 

2.  The  initial  diameter  of  the  vortex  core  was 
about  10  percent  of  the  blade  chord  and  more  than 
doubled  its  size  after  one  revolution  of  the  rotor. 
The  rate  of  growth  followed  a  square-root  function 
of  wake  age. 

3.  The  boundary  layer  wake  that  was  shed  along 
the  span  of  the  rotor  blade  disappeared  after  ip  = 
30°;  however,  two  concentrated  zones  of  oppositely 
directed  axial  flow  were  observed  near  the  center  of 
the  vortex  throughout  the  first  rotor  revolution. 

4.  The  void,  vorticity  (or  centroid  of  vorticity), 
and  the  swirl  centers  occurred  at  distinct  locations, 
but  were  never  separated  by  more  than  26  percent  of 
the  vortex  diameter. 

5.  The  axially  directed  velocity  profile,  relative  to 
the  convection  velocity  of  the  void  center,  was  asym¬ 
metric  because  of  the  influence  of  the  rotor  down- 
wash.  The  vortex  also  appeared  to  have  an  elliptical 
shape  that  rotated  as  the  wake  aged. 


6.  The  peak  vorticity  occurred  near  the  center  of 
the  vortex  during  the  early  development  of  the  trail¬ 
ing  vortex,  and  then  progressively  shifted  toward  the 
core  diameter.  It  is  presently  unknown  to  what  ex¬ 
tent  this  result  is  due  to  weakly  defined  characteris¬ 
tics  within  the  void  region. 

7.  The  magnitude  of  the  radial  distribution  of  the 
circulation  rapidly  increased  until  reaching  a  point 
just  beyond  the  core  radius,  after  which  the  rate  of 
growth  decreased  significantly.  The  core-radius  cir¬ 
culation  increased  slightly  with  wake  age,  but  the 
large-radius  circulation  appeared  to  remain  relatively 
constant. 

8.  The  radial  distributions  of  swirl  velocity  and 
vorticity  exhibit  self-similar  behaviors,  especially 
within  the  core.  According  to  vortex  models  that 
approximate  the  measured  data,  the  core-radius  cir¬ 
culation  was  about  79  percent  of  the  large-radius  cir¬ 
culation,  and  the  large-radius  circulation  was  about 
67  percent  of  the  maximum  bound  circulation  that 
was  calculated  for  the  rotor  blade.  On  average,  about 
53  percent  of  the  maximum  bound  circulation  resides 
within  the  vortex  core  during  the  first  revolution  of 
the  rotor. 
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Table  2:  Locations  related  to  vortex 


Wake  age 
(deg) 

Vorticity  center 
r  (mm)  z  (mm) 

Swirl  center 
r  (mm)  z  (mm) 

Void  center 
r  (mm)  z  (mm) 

Core  radius 
(mm) 

Void  radius 
(mm) 

0 

945.78 

2.89 

946.11 

1.93 

946.00 

3.00 

5.00 

0.71 

15 

934.78 

-1.01 

934.25 

-1.31 

934.70 

-0.80 

5.25 

2.21 

30 

922.28 

-7.88 

922.99 

-8.66 

922.10 

-8.80 

5.00 

2.85 

45 

909.53 

-16.08 

909.51 

-14.88 

910.00 

-16.00 

6.75 

3.42 

60 

899.23 

-26.00 

899.35 

-25.80 

899.30 

-26.70 

6.75 

3.74 

75 

888.40 

-35.37 

888.98 

-34.31 

888.40 

-35.20 

6.75 

4.42 

90 

877.85 

-44.60 

878.04 

-45.51 

878.20 

-45.20 

7.50 

5.04 

105 

868.99 

-54.56 

867.06 

-54.17 

869.00 

-55.20 

7.75 

5.44 

120 

858.57 

-64.35 

858.66 

-63.02 

858.50 

-64.50 

8.50 

5.55 

135 

850.34 

-76.28 

851.04 

-73.69 

850.00 

-75.70 

8.50 

5.84 

150 

842.34 

-83.49 

839.95 

-84.16 

842.20 

-85.40 

10.50 

5.94 

180 

829.99 

-112.12 

828.45 

-112.83 

830.60 

-111.80 

8.75 

6.71 

195 

820.94 

-130.12 

820.29 

-130.54 

823.60 

-131.10 

8.75 

7.08 

210 

820.87 

-148.14 

819.50 

-146.08 

821.10 

-149.40 

9.51 

6.80 

225 

814.53 

-165.86 

812.56 

-166.17 

816.60 

-166.00 

8.75 

6.82 

240 

810.73 

-185.65 

808.70 

-183.76 

812.70 

-185.30 

9.99 

6.62 

255 

809.76 

-200.91 

808.47 

-201.46 

810.50 

-202.00 

10.24 

7.49 

270 

810.10 

-218.07 

807.03 

-217.92 

810.70 

-217.60 

9.75 

6.51 

285 

804.78 

-231.66 

801.58 

-231.76 

804.50 

-232.60 

10.00 

6.28 

300 

802.36 

-246.28 

800.52 

-245.89 

803.10 

-247.00 

9.75 

6.57 

315 

802.14 

-262.10 

800.28 

-261.70 

802.40 

-263.10 

9.75 

6.72 

328 

801.18 

-277.29 

803.30 

-276.32 

801.10 

-277.00 

10.00 

7.04 

344 

794.13 

-293.91 

791.73 

-295.04 

794.50 

-293.60 

10.50 

6.76 

375 

795.78 

-312.19 

796.70 

-309.02 

795.50 

-311.90 

11.00 

7.11 

390 

798.17 

-328.52 

797.79 

-325.25 

798.90 

-329.10 

11.25 

7.14 

Table  3:  Vortex  convection  and  shape 


Wake  age 

(deg) 

Vortex  convection 
uc  (m/s)  wc  (m/s) 

Max  swirl  velocity 
(m/s) 

Velocity  peak  separation 
<£/,  (mm)  dv  (mm) 

0 

-3.43 

-3.03 

30.77 

10.77 

9.07 

15 

-3.43 

-3.04 

35.36 

9.58 

11.15 

30 

-3.41 

-3.07 

32.24 

10.93 

11.02 

45 

-3.37 

-3.14 

24.88 

13.64 

12.88 

60 

-3.24 

-3.31 

25.86 

13.51 

14.55 

75 

-3.23 

-3.39 

27.39 

13.68 

12.75 

90 

-3.12 

-3.57 

24.59 

13.65 

14.58 

105 

-2.99 

-3.79 

23.26 

13.67 

14.56 

120 

-2.84 

-4.04 

20.29 

13.49 

18.30 

135 

-2.66 

-4.30 

16.62 

14.91 

18.36 

150 

-2.48 

-4.57 

17.07 

16.26 

18.19 

180 

-2.09 

-5.07 

20.78 

17.61 

14.46 

195 

-1.89 

-5.28 

17.50 

19.03 

18.14 

210 

-1.70 

-5.43 

17.94 

18.04 

17.81 

225 

-1.52 

-5.53 

17.52 

17.99 

15.68 

240 

-1.35 

-5.58 

15.32 

19.37 

17.65 

255 

-1.20 

-5.58 

17.82 

19.38 

19.71 

270 

-1.06 

-5.53 

17.27 

19.02 

18.26 

285 

-0.95 

-5.45 

15.94 

18.99 

16.44 

300 

-0.85 

-5.35 

19.57 

18.97 

18.32 

315 

-0.77 

-5.25 

18.20 

19.04 

18.10 

328 

-0.72 

-5.17 

17.71 

19.07 

18.26 

344 

-0.68 

-5.08 

17.84 

19.02 

18.29 

375 

-0.64 

-4.99 

16.35 

18.96 

21.74 

390 

-0.64 

-4.99 

15.63 

19.11 

21.82 

Table  4:  Vorticity  and  velocity  extrema 


Vorticity  normal  to  Velocity  within  Velocity  normal  to 


Wake  age  measurement  plane  measurement  plane  measurement  plane 


(deg) 

max  (1/s) 

min  (1/s) 

max  (m/s) 

min  (m/s) 

max  (m/s) 

min  (m/s 

0 

20463.3 

-2484.1 

36.32 

1.08 

4.79 

-10.78 

15 

22942.8 

-994.0 

38.76 

1.21 

2.62 

-9.40 

30 

21014.8 

-1290.1 

35.08 

1.54 

1.83 

-9.25 

45 

11374.9 

-1198.1 

26.59 

0.23 

2.33 

-7.61 

60 

11808.0 

-1055.6 

26.67 

0.02 

2.79 

-7.16 

75 

11929.1 

-1043.4 

29.03 

0.14 

2.71 

-6.72 

90 

10505.9 

-918.7 

26.52 

0.05 

4.57 

-6.61 

105 

9416.3 

-1028.0 

25.78 

0.05 

4.49 

-6.56 

120 

7657.5 

-1050.3 

23.93 

0.03 

4.69 

-6.55 

135 

5812.4 

-1172.7 

22.56 

0.20 

5.03 

-6.96 

150 

6990.6 

-1288.0 

25.69 

0.85 

4.55 

-9.30 

180 

8575.8 

-1143.9 

24.23 

1.37 

3.31 

-5.31 

195 

8755.1 

-1814.7 

23.55 

0.72 

3.52 

-6.09 

210 

8154.6 

-900.5 

23.14 

0.61 

3.48 

-5.76 

225 

7746.3 

-1147.3 

22.40 

0.44 

3.26 

-4.50 

240 

6852.9 

-1623.4 

19.61 

0.24 

3.46 

-4.41 

255 

7020.4 

-724.5 

20.62 

0.76 

4.19 

-5.36 

270 

7006.7 

-870.7 

22.13 

0.27 

2.77 

-4.82 

285 

6732.1 

-1318.6 

19.52 

0.06 

2.40 

-4.43 

300 

7062.0 

-706.5 

22.23 

1.37 

2.11 

-5.33 

315 

6382.8 

-773.3 

19.93 

1.20 

2.35 

-3.68 

328 

6532.7 

-788.4 

21.20 

0.90 

2.79 

-4.44 

344 

6841.8 

-859.0 

20.49 

0.38 

1.72 

-9.15 

375 

5851.2 

-492.8 

.20.95 

0.73 

2.63 

-7.30 

390 

5977.2 

-734.8 

21.39 

0.05 

4.32 

-7.26 
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Figure  1:  OLS  model  rotor  blade. 


Figure  2:  Setup  for  PIV  measurements. 
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Figure  4:  Sequence  for  extracting  vortex  characteristics. 
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Figure  5:  Wake  ages  at  which  PIV  data  were  recorded  in  zones  A,  E,  and  F. 
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Figure  7:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ip  =  30°. 
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Figure  8:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ip  =  90°. 
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Figure  9:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  ip  =  195°. 
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Figure  10:  Velocity  and  vorticity  characteristics  with  reference  to  the  center  of  swirl  at  rp  —  375°. 
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Figure  14:  Velocity  and  vorticity  profiles  along  several  cross  sections  of  the  trailing  vortex  at  tp  =  15°  and 
195°. 
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Figure  16:  Vorticity  (w«),  calculated  along  radius  k  relative  to  the  vorticity  center  using  the  inplane  velocity 
field  at  ip  =  15°,  195°,  and  375°. 
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Figure  18.  Vortex  size  and  strength  during  wake  ages  ip  =  0°  to  390°  and  hover  inflow  models. 


